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Summary. — We present a detailed study of the spectral and temporal properties 
of the X-ray flares emission of several GRBs. We select a sample of GRBs which 
X-ray light curve exhibits large amplitude variations with several rebrightenings 
superposed on the underlying three-segment broken powerlaw that is often seen 
in Swift GRBs. We try to understand the origin of these fluctuations giving some 
diagnostic in order to discriminate between refreshed shocks and late internal shocks. 
For some bursts our time-resolved spectral analysis supports the interpretation of 
a long-lived central engine, with rebrightenings consistent with energy injection in 
refreshed shocks as slower shells generated in the central engine prompt phase catch 
up with the afterglow shock at later times 

PACS 98.70.Rz - 7-ray sources; 7-ray bursts. 
PACS 95.85. Nv - X-ray flares. 



1. Introduction 

Most of the Gamma Ray Bursts (GRBs) observed in the pre-Swift era showed a 
smooth power law decay with time in the X-ray and optical afterglow light-curve pQ. 
This behavior was found to be consistent with the basic predictions of the fireball model, 
where the afterglow flux is produced as a relativistic blast wave propagating into an 
external medium. Under the assumption of a spherical fireball and a uniform medium, 
Sari, Piran & Narayan (1998) showed that the dependence of the flux on frequency 
and time can be represented by several power law segments, F v oc v~$t a , each of which 
applies to a different regime. Before Swift only a few GRBs showing deviations from the 
smooth power law light curve were known, see e.g. the case of GRB 021004 0,0]. 

This simple picture is now changing since the advent of Swift. Erratic X-ray flares 
have been detected by Swift in almost half of its GRBs (0, 0, 0, 0, 0)- While some 
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bursts show one distinct flare after the GRB prompt emission, (e.g. GRB 050406 |10j). 
other events like GRB 050502B and GRB 050713A, show several flares 0, 0, One of 
the main goals of current GRB science is to understand the origin of this newly observed 
light curve behavior. Since flares are likely to trace the activity of the central engine, 
they can help us gain insights on the physical processes governing the early post-GRB 
phases. Note that flares have been observed in the light curves of both long and short 
GRBs. As a first step toward the understanding of different flares characteristics it is 
fundamental to carry out a systematic investigation of the morphological and timing 
properties of the flares observed by Swift/XRT (Chincarini et al. 2006J2HJ ) ■ 

In some cases we can already rule out some models that have been proposed to explain 
the flare phenomenon. The presence of several flares in one GRB strongly disfavors 
scenarios that can only account for one flare, e.g. the self synchrotron Compton emission 
from reverse shock Density bumps |12| cannot be responsible of large amplitude 

rebrightenings (e.g. GRB 050502B 0) . In this paper we carry out a detailed spectral 
and timing analysis of flares in several GRBs with the goal of constraining the physical 
mechanism responsible for its production. In particular we concentrate on two models 
proposed for the flares, the late internal shock (IS) model and the refreshed shocks (RS) 
model. 

In the internal shock scenario the flares are due to the reactivation of the GRB central 
engine and they are produced from an "internal" dissipation radius within the external 
blast wave front In this case a mechanism to reactivate the engine is required. For 
the long bursts, King et al. (2005) ^3] proposed a model in which the flares could be 
produced from the fragmentation of the collapsing stellar core in a modified hypernova 
scenario. For the case of short GRBs, MacFadyen et al. (2005) suggested that the 
flares could be the result of the interaction between the GRB outflow and a non-stellar 
companion. More recently, Perna et al. (2006) [Fg\ have analyzed the observational 
properties of flares in both long and short bursts, and suggested a common scenario in 
which flares are powered by the late-time accretion of fragments of material produced 
in the gravitationally unstable outer parts of an hyper-accreting accretion disk. Other 
mechanisms that could produce flares involve a magnetic origin |17j . |18j . 

In the refreshed shock scenario the flares are due to late time energy injections 
into the main afterglow shock by slow moving shells ejected from the central engine 
during the prompt phase. In this case we do not need a reactivation of the engine. This 
model can interpret several flares features [301, [21] seen in some GRBs. 

Both models have observational consequences that can be tested. In this paper we 
propose and exploit some diagnostics between the two models. 

2. The sample and the data analysis 

We consider all the GRBs detected by Swift between the Swift launch and 2006 
January 31. We select all those that showed prominent flares in their X-ray lightcurves. 
In table 1 we list all the GRBs that were selected for the analysis. A more extensive 
analysis that will be presented includes a larger sample of Swift GRBs The X-ray 
data were reduced using the Swift Software (v. 2.0) and in particular the XRT software 
developed at the ASDC and HEASARC [23] ( x ). 

Important information on the mechanism producing the flares can be obtained from 
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the comparison between the observed variability timescale AT and the time T, at which 
the flare is observed. In the refreshed shock scenario the early afterglow has a variability 
timescale of the order of the time since the explosion, AT ~ T JHj; while in the late 
internal shock model AT << T. We calculated for all the flares of the sample the ratio 
between the duration of the flare, AT, and the time at which the flare peaks T pea k and 
report these values in Table 1. We also report the time at which the rising behavior of 
the light curve begins T r i so - In our analysis AT is defined as the difference between the 
times at which the light curve shows the end of the falling and the beginning of the rising 
behaviors. All times refer to the BAT trigger. As we can see from this table, within the 
uncertainties, all flares have AT/T pea k ~ 0(1), in agreement with the refreshed shock 
scenario. However some of these flares have other characteristics that can be better 
explained in the internal shock scenario. 

The strategy for our temporal fitting procedure is the following. We consider the data 
in the range between T r j SC and T pca k to compute the slope, a r , of the rising part of the 
light curve. The start time is set to T r i se . For the decay, we use the data in the range 
between T poa k and the end of the decaying behavior to estimate the power law decay 
index ad- The start time has been set to T poa k. We note that the power-law rise and 
decay indices depend strongly on the initial counting time, To. In the internal shock 
scenario, the re-starting of the central engine re-sets the starting time. This is why some 
authors [21] use Tq ~ T r j SC to estimate ad and thus find much steeper temporal slopes 
(we call these slopes aa,is)' We think that the right choice for the start time is T pea k. 
This is because a start time earlier than the beginning of the decaying behavior causes 
an unphysical steepening of the slope, due only to mathematical effects. In any case, 
since the To issue is an open problem, we also discuss what happens to the decay slope if 
the start time is set to T r i se . We therefore report both ad and ad, is for the bursts were 
the internal shocks better interpret the flare phenomenology. 

2T. Time-resolved spectroscopy. - The analysis of the window timing (WT) light 
curves of the flares reveal complex spectral variations. To investigate the nature of these 
variations we performed a time resolved spectral analysis. We first fitted the spectra in 
the rise and decay parts with a simple power law combined with photoelectric absorption. 
We indicate with T = (3 + 1 the photon index. The results of our fits are shown in Table 

2. Large spectral variations are evident, in particular between the spectra corresponding 
to the rise and fall of several flares. 

Motivated by the synchrotron emission model 0, we then fitted the spectra using a 
broken power law model. We kept the low energy (before the cooling frequency) spectral 
power law index fixed at Pi = T± — 1 = 0.5, as expected from standard synchrotron 
models and left the high energy power law index /?2 = — 1 as a free parameter. The 
results are shown in Table 2. In some cases we were not able to fit the spectrum with a 
broken power law, as evident in Table 3 (we write "na" in this case). 

3. Refreshed or late internal shocks? 

In this section we try to distinguish between the two (IS and RS) models using the 
temporal and spectral properties derived for each flare in the previous section. As we 
have emphasized, the parameter that helps in this analysis is the flare duration, AT. 

3T. Late-time "energy injection" . - If AT ~ T rcbrightenings and bumps in the 
afterglow light curve may be due to "energy injections" at later times that could be 
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Table I. - values of z, AT, T for each flare 



GRB 


z 


T ■ 

1 rise 


-*peak 


AT 


AT/Tpeak 


050406 




132 


212 


160 


0.75 


050502B 




468 


848 


700 


0.83 


050713A 




97 


107 


50 


0.47 


050730(1) 


3.97 


202 


232 


100 


0.43 


050730(2) 


3.97 


302 


432 


300 


0.69 


050730(3) 


3.97 


602 


672 


340 


0.51 


050822 




211 


231 


70 


0.30 


050904 


6.29 


368 


458 


200 


0.44 


050908 




268 


388 


500 


1.29 


060111A 




204 


274 


260 


0.95 


060124(1) 




311 


571 


330 


0.58 


060124(2) 




641 


701 


230 


0.33 



produced by slower shells that catch up with the afterglow shock at later times. In this 
case we have some predictions on how the temporal and the spectral energy slope should 
be related. We consider the standard afterglow model in the adiabatic evolution 0. 
The analysis can be easily extended to the radiative evolution and we will do this in our 
future paper |22|. We can distinguish four main cases that we show in Fig. 1: 

• (i) The spectral break frequencies Vb, (where we have the change in the spectral 
slope) are within the 0.5-10 keV band covered by the XRT observation, both during 
the rise (vb,r) and the decay (vb,d) phases. In this case ad — —1/4 for v < Vb.d 
and ad — — 3p/4 + 1/2 for v > Vb,d- Since the spectral index in the decay, 0d, is 
expected to be f3d — —p/2 after the break (as expected in the synchrotron model), a 
relation between the spectral and temporal slope as, ad = 3/2/3^ + 1/2, is expected. 
Indeed, for a completely adiabatic evolution, one has Vb oc E 1 / 2 (where E is the 
burst energy), while for a fully radiative evolution Vb oc E 4 / 7 , both implying that 
the flux rebrightening F2/F1 scales with the frequency change factor as iyb,rl v b,d) ■ 
As an example of this class we consider the GRB 050713A [20] and in particular 
the first flare detected by XRT. In this case ad follows the model prediction as 
shown in Fig. 2; it goes from ad = —0.25 for v < Vb.d, to ad = —1 for v > Vb.d- The 
predicted rebrightening is F2/F1 ~ (vb,r/vb,d) 2 ~ 10, similar to the observed value. 

• (ii) The synchrotron frequency Vb, r is outside of the band during the rise Vb,r > 10 
keV, as implied by the single power-law fit to the data and Vb,d is within the 0.5-10 
keV band covered by the XRT observation. In this case ad = —1/4 for v < Vb.d 
and a d = 3/2/3 d + 1/2 for v > v b ,d and F 2 /F 1 > (10/^, d ) 2 . 

• (iii) The synchrotron frequencies Vb, r is outside of the band during the rise i>b.r > 10 
keV, and shifts below the lower limit of the energy band Vb,d < 0.5 keV during the 
decay. In this case a d = 3/2/3 d + 1/2 and F 2 /F 1 > (10/0. 5) 2 . 

• (iv) The synchrotron frequencies fb,r is outside of the band during the rise and 
decay phase Vb,r < 0.5 keV. In this case ad = 3/2/3^ + 1/2, but we do not have any 
constraint on FijF\. As an example of this class we consider GRB 050730 j2J an d 
in particular its second flare shown in Fig. 3. 
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Case i 



Raise 



Decay 



-1/2 






\-p/2 



0.5keV v br iokeV 0.5keV v bd lOkeV 



Case ii 

Raise 



Decay 



^1/2 






\ -p/2 



0.5keV 10kcVV br 0.5keV v bd lOkeV 



Case iii 
Raise 



Decay 



-p/2 



0.5keV lOkeV V br V M 0.5keV lOkeV 



Case iv 



Raise Decay 



-p/2 



-p/2 



v br 0.5keV lOkeV v b il5keV lOkeV 



Fig. 1. - The four cases in the refreshed shock model explained in sec. 3.1 



From the best fit spectral energy index during the decay, /3d = 0.70 ± 0.07, we 
find ay = —0.55 which is consistent with the temporal decay index found from 
the temporal fit, aa = —0.67 ± 0.16 as we can see in Table 2. Also the flare of 
GRB050904 satisfies these relations. 

3'2. Late internal shocks. - (AT << T ) 

A diagnostic to check if the re-brightenings are due to late internal shocks has been 
recently proposed by Liang et al. 2006 24 . In the internal-origin scenario for X-ray 
flares the decay emission of the flaring episodes should be dominated by high-latitude 
emission with the temporal decay index related to the decay spectral index as ad = /3d + 2 
f |25|h Liang et al. 2006 use the time integrated spectrum; however in order to test the 
curvature effect, the analysis should be done using the spectrum in the decay phase of 
the flare. We have checked this interpretation for the first three bright X-ray flares 
observed in GRB050730 afterglow shown in Fig. 3 and we plan to extend the analysis to 
all flares in our sample [221 • We used the spectral indices measured during the decay 
and found best fit zero times [21] t ,i = 77±H s, i ,2 = 234±gg s and t ,3 = ^4til 2 s : 
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Fig. 2. - Power-law decay indices of the GRB050713A second flare light curves in the 0.3-0.8, 
0.8-1.4, 1.4-2.2, 2.2-4, and 4-10 keV bands, starting from 108 s after to, as a function of the 
energy. The prediction of the synchrotron model is represented by the dashed lines, considering 
a broken power-law for two different values of p. 



respectively. Comparing these values with the rise times given in Table 1, we can see that 
for all three episodes the zero time values are not located during the rising portion of the 
corresponding flare, clearly indicating that the observed decay slopes are not consistent 
with being due to high-latitude emission, as predicted by the late internal shock scenario. 

Nonetheless there are flares that can be better explained with the internal shock model 
even if AT - T. These are for example GRB050406 and GRB050502B where the presence 
of an underlying continuum in the X-ray lightcurve, consistent with the same slope before 
and after the flare argues against the external shock models in the adiabatic regime, since 
no trace of an energy injection can be found. However it is possible that for these flares 
the radiative regime applies |22| . If we consider To = T r i se we estimate ay, is = — 3.2±0.7, 
a dIS = -4.3±0.1 and a d ,rs = -1.4±0.2 for GRB050406, GRB050502B and GRB050822 
(that is the only one that has AT << T). We can compare these values with the ones 
expected by the curvature effect. For GRB050406 we find ad, is = — {fid + 2) = —3.2, in 
agreement with the observed value, whereas for GRB050502B and GRB050822 we find 
oy.is — —3.3 and ad = —3.8, which are different from the observed values. 

4. Discussion 

It is by now apparent that flares in the light curves of GRBs decay have different 
characteristics. While it may prove difficult to explain their complex phenomenology 
within the framework of a single model, we can nevertheless constrain some of the models 
by carefully exploiting the present data. 

We have presented a timing and spectral analysis of several flares detected by Swift. 
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Fig. 3. - Swift XRT 0.3-10 keV light curve of the GRB 050730 X-ray afterglow during the 
first orbit. The solid line is the best fit model to the data obtained considering a linear rise 
exponential decay for the three flares (see Sect. ?? of Perri et al. 2006). The dashed vertical 
lines delimit the seven time intervals considered for the spectral analysis. Data are binned to 8 
seconds resolution and errors are at the lcr level. 



Since the main energy output of the GRB engine rapidly decays with time |26|. the 
energy injections are likely due cither to slower shells ejected during the prompt phase 
of the GRB engine and catching up at later times (refreshed shocks), or to later energy 
production episodes in the central engine, due, e.g., to disk fragments that accrete at later 
times (late internal shocks). We have proposed a method to distinguish between these 
two different models based on the estimate of the flare duration and on the comparison 
between the spectral and temporal slopes. 

We found that the flares of GRB050713A, GRB050730 and GRB050904 can be ex- 
plained by the refreshed shock model. We saw that the flares in GRB050730 do not 
satisfy the relation implied by the curvature effect between the temporal and spectral 
slope, indicating that internal shocks may be less likely in this case. However there are 
other bursts (like GRB 050406, GRB 050502B and GRB050822) were late internal shocks 
may interpret better the flares characteristics. We plan to extend our analysis to a larger 
sample of Swift GRBs. 

Flaring activity is also observed in the optical afterglow of some GRBs (like GRB 
050730 see Perri et al. 2006). Detection of optical rebrightenings in the light curve 
would also help in testing the models. Simultaneous re- brightening events in the X-ray 
and optical bands, like the ones seen in GRB 050730 can be easily explained in the 
refreshed shock model (e.g. |2"?]). 
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Table II. - Best fit values for the temporal and spectral power slopes. We also 
indicate which models best interprets the flares phenomenology, internal shocks 
(IS) and refreshed shocks (RS). For the RS model we also report to which case of 
the ones discussed in sec 3.1 the GRB belongs 



GRB 


a 




r 


xl 


model 






10 22 cm~ 2 






050406(r) 


1.1 ±0.4 


< 0.09 


2.45±0.71 


0.6 




050406(d) 


-1.05 ±0.3 


< 0.06 


2.23±0.55 


1.42 


IS 


050502B(r) 


1.41 ±0.03 


0.15±0.01 


2.27±0.03 


1.02 




050502B(d) 


-0.88 


0.08±0.01 


2.33±0.04 


1.2 


IS 


050713A(lr) 


0.67 ±0.05 


0.53 ±0.09 


1.60 ± 0.15 


1.11 




050713A(ld) 


-0.25 to -1 


0.62 ± 0.05 


2.60 ±0.15 


1.27 


RS case i 


050730(lr) 


0.37 ±0.12 




1.29 ± 0.16 


0.86 




050730(ld) 


-0.41 ±0.06 




1.82 ±0.12 


1.3 




050730(2r) 


0.92 ±0.10 


o i +0.8 
z - L -0.7 


1.71 ± 0.08 


1.1 




050730(2d) 


-0.67 ±0.16 


0.9 ±0.5 


1.70 ±0.12 


1.1 


RS case iv 


050730(3r) 


0.51 ±0.12 


u -'-0.7 


1.77 ± 0.12 


0.87 




050730(3d) 


-0.27 ±0.08 


1.0 ±0.6 


2.01 ±0.1 


0.81 




050822(r) 


0.45 ±0.15 


0.24±0.05 


2.8±0.2 


0.8 




050822(d) 


-0.56 ± 0.08 


0.19±0.04 


2.8±0.2 


0.9 


IS 


050904(r) 


0.25 ±0.06 


0.14±0.02 


1.76±0.06 


1.22 




050904(d) 


-0.64± 0.05 


0.07±0.03 


1.75±0.08 


0.8 


RS case iv 


060111A(r) 


0.000 


0.14±0.02 


1.76±0.06 


1.22 




060111A(d) 


0.48± 0.03 


0.07±0.03 


1.75±0.08 


0.8 


RS 
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GRB 


10 22 cm~ 2 


r 2 


E 6 


xi 


050406(r) 


na 


na 


na 


na 


050406(d) 


na 


na 


na 


na 




fl 1 94- fl f)9 


9 9^4- fl 04 


n 71 + n ns 


u. yu 


050502B(d) 
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2.29 ± 0.05 


0.72± 0.08 
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050713A(lr) 


0.50 ± 0.09 


1 - yo -0.35 


J -1.75 
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050713A(ld) 
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050730(lr) 


1 6+ 2 ' 5 
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2 n +11 
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1 l +0 - 7 
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na 
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na 


na 


050730(2d) 


na 


na 


na 


na 


050730(3r) 


na 


na 


na 


na 
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na 


na 
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na 


na 
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0.8± 0.7 


0.8 
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